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A path integral over trajectories of 2n fluid particles is identified with a 2n-th order correlation 
function of a passive scalar convected by d-dimensional short-correlated multi-scale incompressible 
random velocity flow. Strong intermittency of the scalar is described by means of an instanton 
calculus (saddle point + fluctuations about it) in the path integral at n 3> d. Anomalous scaling 
exponent of the 2n-th scalar's structural function is found analytically. 
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pH i 1. INTRODUCTION 

^ : 

'"^ i The problem of scaling behavior in Kraichnan's model of a white-advected passive scalar [Q, attracts a great deal of 
; attention |,U4||,|,g. In the wide range of scales, called the convective interval, structural functions of the Lagrangian 
i tracer 9, passively advected by d-dimensional short-correlated in time multiscale incompressible flow, possess a scaling 
I ' behavior. The anomalous scaling exponent (^2n of the 2n-th order structural function, {{9i — 6*2)^") oc r^^", has been 
J> . calculated in the following cases: i) large space dimensionality ^ (2 — C2)"-, for n = 2 in [Q, and generally for 
T— I ' all allowed n in 0, (2n 2n{n — 1)(2 — C,2)/d'i ii) almost smooth scalar field 2 — C2^1,c?>2 for n = 2 in § and 
generally for n(2 — (,2) ^ d, d > n va. C,2n 2n{n — 1)(2 — C2)/(rf + 2). The perturbation methods yield the scaling 
exponents in the limits where the respective bare approximations are strictly Gaussian and the anomalous corrections 
are small. 



o 
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' Instanton (steepest descent) formalism after perturbation expansion is the second quantitative method that could 
0^ ' be applied to a general statistical problem. The method works when a large parameter makes some very special rare 
' configuration to have an exponentially large weight. Such a large parameter may be a high order n of correlation 
^.^i function ((^") of fiuctuating field (p. The bare instanton approximation is obviously strongly Non-Gaussian. The 
''O ' idea was originally introduced and successfully applied in field theory almost twenty years ago, but introduced 
Q , to the turbulence theory only very recently. An instanton calculus in a Lagrangian path integral was used to find 
' an exponential tail of the scalar's probability distribution function (refiected intermittent, non-Gaussian behavior of 
,.Ch . higher moments) in the case C2 = of linear velocity profile |^ (later on it was shown that the limit turns out to be 
J-^ solvable exactly ) . A general method for finding the non-Gaussian tails of pdf for solutions of a stochastic 

J> differential equation, such as the convection equation for a passive scalar, random driven Navier-Stokes etc., was 
formulated in . The initial idea of the method is to look for a saddle point configuration in the path integral for 
the generating functional introduced in . The extremum of the effective action is given by a coupled field- force 

' configuration (instanton), varying in space and time. The method was applied recently to Burgers' turbulence Jl7| , p^ . 
Generally, it is very difficult to solve the coupled (field- force) instanton equations. 

In the present paper we generalize the idea of |^ for the case of a nonlinear velocity profile (C2 > 0) . C2n is calculated 
for n being the largest number in the problem. The method is based on a very special feature of the problem [p^p[ : 
there exists a closed differential equation connecting 2n-th and (2ri — 2)-th simultaneous correlation functions of the 
scalar. The 2n-th correlator is expressed via the convolution of the resolvent of the eddy-diffusivity operator with a 
source function constructed from the (2n — 2)-th correlation function. To prepare a path- integral for the instanton 
calculus, we perform an explicit map of the original problem of calculation of the 2n-th order correlation function 
to the problem of calculation of a matrix element in an auxiliary 2n-particle Quantum mechanics. The resolvent of 
the eddy-diffusivity operator is expressed in the method via the path integral over trajectories of 2n fluid particles 
moving from an initial geometry (at which we are aimed to describe scalar's correlations) with a characteristic scale 
r to a final large-scale (~ L) geometry. The tensor of eddy diffusivity plays the role of tensor of inverse mass for the 
particles from the associated Quantum mechanics. The tensor depends explicitly on relative distances between the 
particles. 

It is the large number of particles that makes the auxiliary Quantum mechanics almost "classical" ( "semi-classical" ) . 
A "classical" 2n-particle configuration is the desirable rare event that describes both the intermittency of 2n-th 
moments of scalar differences and intermittency of the n-th moment of the dissipation field e = k(V6')^ (it is proven 
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in 1^ that they are related to each other, if scale-invariance of the structural function and of the correlator of the 
dissipation field is valid: (e") ^ [L/r^]'^^", A2„ = 71(2 — (271', k and are the diffusion coefficient and scale, 
respectively). Calculation of the "classical" (saddle-point) contribution into (e") gives the scale- invariant answer: 
^2n ~^ '^C2 at n — > 00. The "classical" anomaly shows the highest level of intermittency possible: = 0- Therewith 
exist a wide set of "classical" trajectories (realizing themselves separately, for different initial displacements of the 
points and different forms of the scalar's source) responsible for the "classical" answer. To extract an optimal trajectory 
(that gives the lowest possible contribution of fluctuations) from the set of the "classical" ones and thus to get a finite 
asymptotic for C,2m we must account for the fluctuations about the saddle-points. 

It is shown that, at < C2 < 1, the optimal trajectory is defined as a relative dispersion of two groups (drops) of 
particles: there is one distance (separation between the drops) being stretched while all the other distances (sizes 
of the drops) are being contracted dynamically. They are Gaussian fluctuations about the optimal trajectory that 
should give the true value of (2n- Accounting for the relative longitudinal (along the "classical" stretching direction) 
Gaussian fluctuations of the drops, gives the dominant contribution into the n-independent asymptotics for (^2n at 
n 00. The exponent is finite, it grows linearly with d and decreases monotonically with an increase in C2. The 
finite limit for (211 at C2 — > along with the nonanomalous answer C,2n = for the strictly "logarithmic" limit C2 = 
show together a discontinuity of (^2n at C2 = 0. There exists a simple physical picture that explains the origin 
of this discontinuity. In the first case of a linear velocity profile, distances between all the fluid particles are stretched 
homotopically: there is no way for two groups of particles to diverge from each other and to keep the inner group 
distances contracted (or even intact) simultaneously. Vice versa, in the case of finite C2 (yet (2 should be smaller than 
unity), the two-point trajectory, with the sizes of the drops of particles being contracted dynamically whereas the 
distance between the drops being increased, is allowed. 

Most of Non-Gaussian fluctuations about the saddle-point conflguration can be dropped in comparison with the 
Gaussian ones if n 3> Pe, d (we should only worry about the explicit calculation of the Non-Gaussian fluctuations 
corresponding to a soft rotation mode). This method is not applicable for Pe being of the order of (the more so as 
being larger than) n. However, making use of an overall observation (concerning the linearity of the problem and the 
scale invariance feature of different terms entered in the correlation functions) one can extend the anomalous result 
(but not the method used for its derivation) to the limit Pe, d too. 

The two-point configuration is not relevant at 1 < ^2 < 2 (repulsion of particles inside of a drop is no longer weak 
to make the configuration stable dynamically). Only trajectories with many (~ n) distances being diverged should 
be taken into account. However, calculation of fluctuations about such trajectories shows a strong renormalization 
of the saddle-point answer: it is a product of ~ n algebraic terms (each responsible for fluctuation of a distance) 
that makes the contribution of fluctuations competitive with (or even large than) the "classical" value. The resulting 
contribution into (e") is negligible in comparison with the normal scaling term, that always exists. To conclude, the 
instanton calculus is not an appropriate tool in this case. 

The material in this paper is organized as follows. In Section II, after a detailed and formal definition of the problem 
we introduce path integral representation for the 2n-th order correlation function of the passive scalar. We present 
the path integral for (e") too. It completes preparation for delivering an instanton (steepest-descent) formalism for 
calculation of (e") at rt 3> d in the two forthcoming Sections. Saddle-point equations are derived and studied in Section 
III. The contributions of different saddle-points into (s") are calculated (Appendix A) and compared with each other 
in Subsection IIIB. To improve the saddle-point calculations and to extract among the saddle-points an optimal one 
we study Gaussian fiuctuations about the saddle-points in Section IV and Appcndises B,C. The anomalous exponent 
^2n for the optimal saddle-point is calculated there in Section IV. In the concluding Section V we discuss the results 
from the points of view of criteria of their applicability, restrictions imposed, possible generalizations, and comparison 
with other results and methods. 

2. FORMULATION OF THE PROBLEM 

The advection of passive scalar is governed by the following equation 

{dt+vc.\/c.- kA)9 = f, V„fe« = 0, (2.1) 

where f{t;r) is the external source, v(t;r) is the advecting d-dimensional velocity and k is the diffusion coefficient. 
/(<; r) and v(t; r) are independent random functions oit and r, both Gaussian and (5-correlated in time. The source is 
spatially correlated on a scale of the pumping L, i.e. the pair correlation function {f{ti; ri)/(t2; r2)) — ^(ti — t2)x(''i2) 
as a function of its argument decays on the scale L. The value of x(0) = P is the production rate of 9^. The eddy 
diffusivity tensor IC"'^ , that describes the Gaussian velocity correlations, 
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{v"ih;r,y{t2; ra)) = Sih - h) [VoS'^" - /C"^(ri - v^)] , (2.2) 

JC^^ir) = [{d+l- -f)S°'Pr^ - (2 - j)r"rf^] , (2.3) 

depends on two parameters: D that defines the level of turbulence and 7, < 7 < 2, that measures a degree of 
nonsmoothness of the velocity field. 



Averaging (2.1) over the statistics of u{t; r) and f{t; r), one gets the closed equation for the simultaneous correlation 



functions of the scalar Fi.2n = (6'(ri) • • • 6'(r„)) 

— ^2nFi...2n = Xl---2n, (2-4) 

Xi-2n = Xi2-F3--2n + pcrmutatious, (2.5) 

n n 

£„ = -J2 l^^^{r, - r,)V?vf + A,;. (2.6) 

i^j i 

The dependence of the source function X2n{Tij ~ L) on i at r^j i is estimated as ~ upscales from L the 

function decays algebraically fast. It is the major information about X2n required for further consideration. 



The equation (2.5) for the pair correlation function [n = 1) was solved explicitly The pair correlator in the 
convective interval, ^ ri2 <^ L, where r'^~'^ = 2(2 — ^)K/[D{d — 1)], gets the following form 

Thus, the pair structural function is shown to have a simple scaling behavior in the convective interval, {{9i — 02^) ~ 
Pr"^ jD, C,2 = 7, to provide the constancy of the flux of 9"^ there. The scaling exponent of C is —7, the function xi2 
does not depend on ri2 deep inside the convective interval so that its exponent is 0, the solution of ( |2.7D thus may 
be presented in the form Fforc + Z, where we separated the so-called "forced" part of the solution (with the scaling 
exponent 7) from the zero mode (that is constant in this case) . It is the forced part that contributes the second order 
structural function. The separation for "forced" terms and zero modes is valid for higher order correlation functions 
as well. It has been recognized independently by the authors of |^,^,|l9| that there are zero modes Z that may provide 
for an anomalous scaling. A zero mode, possessing the slowest downscale decrease among the ones built on 2n points 
(that is not reduced to a sum of zero modes each built on a less number of points), gives the major contribution 
into the 2n-th order structural function of the scalar, for n > 1 Scaling of such a zero mode should grow with n 
to provide the convexity of C,2n as a function of n (it is an immediate consequence of the Holder inequality, see for 
example [|o|). There are two Gaussian limits where there is no anomaly and it is easy to make a classification of zero 
modes of operator C there: a) limit of large space dimensionality, d = 00; b) so-called "diffusive" limit of the smooth 
scalar field, 7 = 2~, (to be precise it was done even in a more restrictive case, when D/(2 — 7) is finite). It was a recent 
breakthrough in the analytic theory of turbulence, when the anomalous exponent C,2n was calculated perturbatively 
in the leading non-Gaussian order in the respective small parameters: 1/d in and 2 — 7 in One emphasizes 

that both the perturbative techniques do not work for sufficiently large moments 2n, when the anomalous corrections 
are of the order of the normal scaling exponent 717. To deal with C,2n for the largest moments, we shall deliver a 
nonperturbative instanton technique. 



The basic equation (2.4) can be rewritten in the following evolution form (one step back from the derivation of 
(O) presented in 0) 



l---2n 



POO 




/ dT exp 


H 


/o 





Xl---2ri 



^ dT Jl[dIl,TZ{T;r,,-R,}xi-2n{'R,}, (2.8) 

here and everywhere below summation over the repeated particle and dimensional indexes will be assumed. 
TZ{T; r,;, R,} is the resolvent of the operator £2n, 

(a* - £2„{r}) n{t; r„ RJ = 5(t) [] J(r, - R,). (2.9) 
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Considering the differential operator under the exponent from the first hne of ( |2.§| ) as a Hamiltonian of a 2n-particle 
Quantum mechanics we can rewrite the resolvent in the Hamiltonian form of the standard Feynman-Kac path integral 



7e(r;r,R) 



p.(T)=Ri 2ti 

l[Vp^{t)VMt) exp [-S{p{ty,p{t)}] 

P,(0)=r, 



pfpf 



where [/C] is defined as: 



[IC] 



af3 



(2.10) 
(2.11) 

(2.12) 



Retarded regularization of the "mass" ([/C]-term) in the action is considered, that means the following discretization 
procedure: tk = fee, e = T/M, k = 0,---,Af, Vp^it) = Ditti^ ^Pi(^fe)' Pii^o) = r*, PiitAi) = Rj, 'Dpiit) = 
UkLi dpi{tk), M oo, 



M-l 



s = 



-pfih+,){pf{t,+,)-pfitk)) 



(2.13) 



fe=0 



where the path integral for the associated Quantum mechanics could be understood as explaining a random (Brownian) 
motion of 2n particles possessing a very special dependence of the tensor of inverse masses [/C] on displacements of all 
the particles. The resolvent represents the probability for the 2n fluid particles to diffuse from the initial geometry to 
the final one for time T. Notice that another 2n-particle representation [ pT| was used to analyze the pumping-free 
(decaying turbulence) two-dimensional case of a linear (7 — 0) anisotropic velocity profile. 

The representation (2.1C-2.13) is useless if we aim to calculate the functional integral explicitly: it would reduce 
one to calculation of the resolvent of C, that is already stated as a generally unsolved problem. Our aim is modest, we 
are going to study the higher correlation functions, or many-particles problem [n 3> d) on the language of a "quasi- 
classical" approximation for the associated Quantum mechanics. The large parameter should allow us to evaluate the 
path integral from the integrand of ( p. 10 ) (or its spatial derivatives, see below) in a saddle-point (instanton) manner. 



is not scale invariant. The integrations over R and T in (2.8) give rise to a huge set of zero modes, describing 



not only the 2n-th structural function but all the lowest ones too (for details of the zero mode ideology see ) 
To separate zero mode giving the dominant contribution into the 2n-th structural function, which is subleading in the 
zoo of the zero modes, we suggest an another oblique way of solving the problem. The idea is to use an exact scaling 
relation between the n-th order moment of the dissipation field e = k[V0]^ and 2n-th order structural function of the 
scalar that was proved in Q by means of the ultraviolet fusion rules discovered in Q| : 



if m-9,y-)^r-]{L/n,) 



then 



(2.14) 



if it is known additionally that (e") is scale-invariant. Let us emphasize that the relation (2.14) between structural 
function and respective correlator of e is based crucially on the expected scale-invariance of both the objects. At n, 
considered to be the largest number in the theory, the scale-invariance over r/L does not need to be a priory valid. 
We will construct an instanton for the correlator of the dissipation field itself 



e")-K" lim 



/ dTJY[dR^ne{T■,r„R,}x2n{'R.^} 



Pi(T)=R, 



L ^ 2n n 

7^,{^;r„RJ= / []pp,(i)I?P.(i) [] [P2fe-2(0)P2fc-i(0)] exp [-5{p(i);p(i)}] 

/ON i fe=l 



(2.15) 



(2.16) 



P.(0)=ri 



It is easy to check by means of direct Gaussian integrations that the discretization condition (2.13) reproduces the 
correct gradient structure of the e correlator (the Hamiltonian form of the path integral allows it to be easily checked). 
A kind of pairing of the space indexes in the p ~ p integrand of (2.16) is arbitrary (for example, one could make the 
integrand symmetric with respect to all permutations of all the particles). 

There are two different specifications that we are free to fix in the problem's set. It concerns initial and final 
conditions imposed. The initial condition is defined by the initial geometry. The final condition is defined by the 
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source function Xin- However, scaling exponents do not depend on a concrete form of the x-function in accordance 
with the general zero-mode ideology • One can use the freedom to make an appropriate choice of the initial 

geometry and the source function. 

It is evident that integration over R; in ( 2.15| ) cannot be performed in the saddle-point manner, if the source 



function is, for example, a uniform constant inside the circle R < L: All the values of R satisfied, R'^ ^ DT (the 



rough observation will be improved later on), give comparable weights in the integrand of (2.15). However, one can 
force a particular final geometry R^ ~ L to be preferable, choosing the source function to get a sharp maximum 
about i? ~ L, where R is an average size, say R = \/^^j^R^f]/\2n\. Then one can include the variation over R^ in 
the common variation procedure adding the term — ln[x2n] to the action. The formal trick is justified by the general 
expectation to get the dominant contribution into the e correlator from a zero mode of operator C It is the universal 
scaling of a zero mode that defines universal (independent on a concrete shape of X2n) scaling of the e correlator. 



Thus, we are going to raise both the p — p and source terms from the integrand of (2.16) into the exponent to 
variate hereafter the effective action 

n 

'^eff = ^ X! [KP2fc-2P2fc-i] ~ ln[x2n], (2-17) 

k=l 

over all the allowed trajectories (over pi{t),pi{t) for all the t from < t < T) in the next Section. 



3. SADDLE-POINT APPROXIMATION 



An instanton is defined from the saddle-point approximation - extremum of the effective action (2.17) with respect 
to fluctuating coordinates Pi{t) and momenta Pi{t) of all the 2n particles: 



dpf{T) ' 



(3.1) 
(3.2) 



where summations over the particle j index and repeated spatial indexes are supposed; j and j* are indexes of 
conjugated particles from a pair (say 1 and 2 or 2n — 1 and 2n); 



d 



D 



„Q ^rj 



The discrete variant of the instanton equations is 

e[IC]tf{h)p^{h+i) + P?{tk) - Pf ih+i) = 0, 

.a , P"(e) 



Pj(e)p,-(e) 

91n[x2 



p?ih), 



M 



dpfitu)' 



(3.3) 

(3.4) 
(3.5) 

(3.6) 
(3.7) 



where k temporal index in ( |3.4[ - |3.5[ ) is running from 1 to M — 1; ti = e = 0+. The equations (3^,3/7), appearing from 
dSpf^/dpf{ ti) — 0, a nd d Spff/dpf{tM) — respectively, explain the rule of parameterization of the (5-functions from 
the rhs of ( |3.2| ) and ( ^3.1]) . To study the saddle point trajectories at the fixed initial geometry r.^ and a fixed form of 
the source function X2n one should solve the following classical equations of motion 



0, 



in the boundary conditions 



ft(0)=r^, ft(T) = R„ 



(3.8) 
(3.9) 

(3.10) 
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where = Pj(e) is related to and the initial momentum pi(e 0) via ( |3.6|); = Piitn) depends on the final 
momentum Pi(tM) via (3.7). Therefore, the problem is reduced to resolving ( |3 .81 ,3.9) with the boundary condition 



(3.10) fixed and with the constraints (3.6 |3.7| ) imposed afterwards. 

The classical Hamiltonian equations of motion (3.8 3.9) possess the standard set of integrals of motion. First of all, 
due to the independence of the Lagrangian (the integrand part of the action S) on time, the energy (that coincides 
with the Lagrangian) is the conserved quantity 



E : 



const. 



(3.11) 



Second, due to the invariance of the action with respect to a uniform shift of all the particles, the momentum 
P" = X^ jPf I is conserved too. Choosing the appropriate coordinate frame, associated with the center of inertia of the 
system of particles, one forces the full momentum P to be equal to zero. Third, due to the invariance of the action 
with respect to uniform rotation (say, around the center of inertia, where = 0) the angular momentum 



(3.12) 



which is a d — 2-dimensional antisymmetric tensor (£"i '"£i is the d-dimcnsional absolutely antisymmetric tensor), is 
the last globally conserved quantity. 



A. Semi-classical analysis for correlators 



Let us consider a particular instanton solution describing dynamical dispersion of particles from a geometry with 
Pij r' at the initial moment of time 0^, to a final geometry of a common-type; with at least one of the distances 
Pij{T) being of the order of R ^ L. The major saddle-point contribution into the s correlator is 



lim 



poo 

Jo 



(3.13) 



where ^^j-j is the effective action Se/f 
any such trajectory is 



(2.17) taken on the classical trajectory pf. The law of temporal evolution of 



p7/2 _ ^'7/2 ^ ^7/2 
i?7/2 



nW 



a{n)TD' 



(3.14) 
(3.15) 
(3.16) 



where 0+ < t < T. The proportionality signs ^ in (|3.l43T6| ) stand to point out that ( 3.14 - 3.16 ) are correct up to 
constant multipliers, d epending on the details of the geometry. Still, there are no extra scale and n dependences in 
the explicit version of ( 3.14 - 3.16 ). T he unive rsal (with respect to the geometry's variation) scaling behavior, p ^ t^/^ , 



follows from equations (any one of) (3.8 3.£ ), by substituting there scaling (over time) ansatz for both p and p fields 



supported by the energy conservation law ( 3.11 ). Here in ( 3.14 - 3.16 ), considering all the distances to be much larger 
than r^i we dropped diffusion for a while. One accounts for diffusion in a special symmetrical case (Appendix Al) 
aiming to show that to get the principal dependence of an ultraviolet divergent quantity on it is enough, generally, 
just to replace all the separations going to by (see also [Q). 

We consider two very symmetrical cases in Appendix A: 1) the uniform expansion of a Sm sphere (2n points 
uniformly distributed on the sphere), the m < d; 2) divergence of two drops, with n_|_ and 2n — particles merged in 
the first and second points (drops) respectively. One can separate all the possible trajectories on two different types 
dependent on how a{n) behaves with n going to oo. Most of trajectories, we will call them "typical", correspond 
to a linear growth of a{n) with n. To specify, the trajectory is "typical" if the volume bounded by a smooth 
{d — l)-dimensional manifold built on the 2n points is not temporarily increased. Relative divergence of the two-point 
geometry (see Appendix A3), the same as expansion of the Sm geometry with m < d — 1, are "typical". An example 
of "non-typical" trajectory is expansion of the Sd sphere (see Appendixes A1,A2 for an explanation of Sm geometry). 
a{n)/n decreases with n going to oo for a "non-typical" trajectory. 
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For a specific kind of source function x (possessing a sharp maximum) chosen, R appears to be L. There 
are two different intervals over the integral time T. First, r' — p(0+) governed by ( |3.6| ) is about r at < T ^ 
L'^/'^r''/'^ /[Da{n)]; For the largest values of the integral time, L^/'^r^/^/[£)a(n)] ^ T, one gets p(0) = 1/r and 
r' ~ [rU^ / {DTa{n))\^/'^'^~'^\ By substituting the saddle-point values, governed by (3.14-3.16) into ( ^.13 ) one 
gets a divergence in the integral at the largest times. The divergence is formal: it should be stabilized by the 
normalization factor, accounting for an algebraic decay of the resolvent with T. We will see below (in Section IV) that 
the algebraic factor in fact come into the game via account for fluctuations to cut the temporal integration in ( 3.13| ) 
at T ~ /[a{n)D]. Thus, for all the values of 7 except some vicinities of the "diffusive" 7 = 2 and "logarithmic" 
7 = limits one gets 



(3.17) 



( 3.17 ) accounts for the principal dependence of (e")'^' on t he Pe clct number only. It is the second interval (of the 
largest values) of T, that give the dominant contributio n in (3.17 ): All the significant dependence on r (or on after 
taking the limit in the rhs of ( t3.13| )) in the integrand of ( |3.13| ) comes from the pg" term via the multiplier r^". Finally, 
accoun ting for the dependence of k, on the diffusion scale results in the anomaly ( 3.17 ). Note, that the saddle-point 
result ( 3.17 ) is generic for all the space dimensions d> 2. 

What is specific about some vicinity of 7^ = 2~ is an expected inapplicability of the saddle point approximation 
there: A growth of 5^^-^ with a growth in n is diminished by decay of the action as 7 goes to 2^. Note, that in 
the naive diffusion lim it, D — 0, (or in the special limit 7 = 2~, D/(2 — j) — const, see [^) a balance between 
different terms in (3.6) differs strongly from the general case: there the r-term can be dropped in comparison with 
r-independent ones, r' turns out to be of the order of the pumping scale L, that results in the absence of anomaly 
(as it should be: the diffusion case is Gaussian). An infinitesimally small deviation of 7 from 2~ in the special limit 
of 1^] makes the po 1/r anomalous solution to be preferable in comparison with the nonanomalous po ^ 1/L one. 
Thus, the 7 — > 2^ limit is indeed very peculiar. The "logarithmic" 7 = limit is very specific too. What is written 
above in the general scheme is valid if {[r'/L]'*/^ — [r'/L]'*)/j <C 1 is satisfied. However, the inequality ceases to be 
true at some very tight vicinity of 7 = 0: The second term from the Ihs of (|3.6|) could be dropped there (the condition 



is opposite to the one which resulted in (3.17)). 

The anomalous answer (3.17) is generic: the scale invariance holds true for the classical trajectories of both "typical" 
and "non- typical" kinds. It accounts for the n-dependent prefactor in the integrand of ( 3.13| ), discriminating between 
the two kinds of instantons. At R and T considered to be fixed, the r-independent term of 5^^^ gets no n-dependence 
in the first case of the "typical" instanton (see Appendix A2,A3 for the two-point instanton and Appendix Al for 
the Sm spherical case with m > d). Vice versa the "non- typical" instanton (it is a spherical case with m — d, 
for example) gets n-dependence from the bare action ( 3.16| ). a{n) goes to zero as n goes to 00 in this case. To 
conclude, the "non-typical" instanton {Sd one) is suppressed in comparison with the "typical" ones. However, we 
cannot distinguish between different "typical" instantons (the Sm instanton with m > d and the two-point instanton) 
on the "classical" level. The suppression of a "non-typical" instanton in d dimensions has a clear physical explanation. 
It follows from the conservation of the volume of a fluid element, prescribed by the incompressibility condition. It 
is a very rare trajectory (that means it has a small weight) that stretches the 2n points forming the Sd sphere and 
conserves simultaneously the volume enveloped by a d — 1 surface built on the 2n points. The surface cannot be 
smooth in this case, it is very fractal. 

Note, that at 1 < 7 < 2 and d = 2 there is not another symmetrical instanton of the type discussed above except 
the S2 one. The two-point instanton that works pretty well at < 7 < 1 turns out to be unstable at 7 > 1: particles 
being initially dropped together into a group try to diverge from each other hereafter. Most probably it is reasonable 
to study another symmetrical instanton with all the particles being elongated into a straight line in this case. Such 
an instanton could be preferable in comparison with the S2 one. We do not yet consider the straight-line instanton 
in the present paper, postponing it for a future study. 

For any solution (already discussed or another) of the auxiliary problems (3.8, ^.S| ,3.1C) one can design such an 
appropriate initial geometry and a particular form of the source function X2n that the trajectory turns out to 
be a unique solution of the full system of the saddle-point equations (3.1,3.2): Fixing and Pi{0) one arrives at a 



unique (due to constraint (3.(;)) initial geometry. Via explicit dynamics and the second constraint (3.7), one finds an 
appropriate form of the source function to make the saddle-point solution self-consistent. 

The anomalous answer (3.17) is scale invariant. In the leading "classical" approximation the anomaly is extraor- 
dinary large: The normal contribution into C271 is fully compensated by the anomalous one, — *■ at n — > 00. 
However, the result gives no possibility for answering the major questions : Is the e correlator scale-invariant at 
n — > 00? And if it is scale-invariant, what is the asymptotic of (^2n at the largest n? To resolve the problems one 
should account for fluctuations. 



7 



It was established in the present Section that the variety of saddle-point solutions (for different and x) gives the 
same scale dependence (3.17). It supports the general statement [^,^ that the dominant contribution into (e") stems 
from a scale-invariant zero mode oi C2n- However, it follows from the same general statement too, that the dominant 
contribution can be lacking for a special kind of source function and initial geometry. It is along this pathways that 
one should optimize the problem with respect to and the x-function, to find the dominant zero mode contribution: 
We must not only find a contribution of fluctuations in (e") but show it is minimal. 



4. ACCOUNT FOR FLUCTUATIONS 



Let us study in the path integral ( 2.16 ), Gaussian fluctuations about yet a nonspecified "classical" trajectory (p'^' 
and p'^' are supposed to be known). The quadratic with respect to fluctuated fields 6p, 6p, a correction to the classical 
action Sf^j is 



1 





-^SpnO)gZ^Sp^{0),+^6pf{T)Cr;6p^iT), 

k 



1 



naf3 _ d"^ ln[x2»] 



dpfdp^J 



(4.1) 

(4.2) 
(4.3) 



where there is no summation over repeated j-indexes in (4.2); the pair of the particles' indexes i and i* describe a 
conjugated pair of particles; JC"^''^ and /C"'^'"'' stand for the first and second spatial derivatives of /C"^ (see explicit 
expression for the first derivative (|]^)). Performing the Gaussian integrations over Sp, Sp one arrives at the following 
expression for the e correlator accounting for the Gaussian fluctuations 



lim 



poo 

/ dTexp[-5,^},]z/' 
Jo 



*0 



i 

Texp 



exp 



dCdt 



Sxi^)), (5x(f) = exp 



2J 



(4.4) 
(4.5) 

(4.6) 
(4.7) 
(4.8) 



Here in (4.4), we use the canonical Quantum mechanical notations for matrix elements. The operator in (4.5) is the 
descendant of the momentum's term from the integrand of (2.16); the diffusive-like state (4.5) is well defined. Texp in 
(4.6) stands for an antichronological ordered exponential. Thus, we came full circle at this stage of the calculations, 
returning back to a problem in the operator representation form (compare the time-ordered exponential from (4.6) 
with the original operator exponent, say from the first line of (2.8)). It follows from (4.6) that '^t can be understood 
as a solution for the following differential equation 



dt + S£ 



(4.9) 



at an initial moment of time '^t{^) = ^'(0;f). 

Let us consider fluctuations about a "typical" saddle-point trajectory with all the distances stretched somehow 
similarly (we will specify the concrete form of the considered instantons later on). Performing rescaling of temporal 
and spatial variables in (4.9) one simplifles it. In the new dimensionless r, variables 
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T 



P 



7/2 



rpi-7/4 



where p{t) is a typically stretched (3.14) classical trajectory, (4.9) gets the following refined form 



dr + C'\ *(t; s,) = S[T-j \n[R/r']j 5x{s^R^-'''^^/T/D). < r < ln[i?/r'], 



'{0-s, = ^DjTR^I^h/r'^~'''^), 



c 



27 



Q/3 

1 



(4.10) 

(4.11) 
(4.12) 
(4.13) 

(4.14) 



where all the new dimensionless matrixes 1C^A,B are time (r) independent. If we exclude divergent degrees of 
freedom (we should worry about uniform rotation of the classical trajectory, that is a soft mode, separately) from 
it becomes a Hamiltonian of a well posed Quantum Mechanics. It is the Quantum Mechanics of I Gaussian oscillators, 
I ^ 2nd. There is thus seen to be a gap in the spectrum of the reduced operator. There are two essential (for present 
consideration) characteristics of the energy spectrum: the value of the gap Ae and the level spacing 5e between 
the ground state and the lowest excited state. Both the energetic characteristics are positive functions of n, d, 7. 
Stretching time r = 7ln[_L/r']/2 is a big parameter due to Pe ^ d. 

There exist two different situations depending on how Se behaves at the largest n. First, tSe is a large parameter 
if Se docs not decay as n grows. Then, it is the evolution of the ground state giving the major contribution 

into ^fy, 



n) 

The multiplier Z^'' is getting smaller with r algebraically 



6x{s.R'-^/^VT/D)j. 



(4.15) 



^c~t — 



exp 



r'V2 



. [//2+Ae]7/[2-7] 

.l) 



det[' 

^2-7/2j. 



(1) 





7Ai3/2 






L 









n-V[2-7] 



(4.16) 



"2, T is considered to be smaller than 



where a typical matrix element of [^~"'^]^ is estimated as [p'^^{t = 0)]^ ' 
L"' /D, and I counts the number of the stretched degrees of freedom. The second possibility is realized if Se is getting 
smaller with n going to 00. Hence it follows that one gets an evolution of a mixed wave packet built from some amount 
of the lowest eigenstates: The wave function of the ground state in (4-15) must be replaced by the wave function 
of the packet. However the multiplier Z^'' is algebraic again and the characteristic size of the packet has the same 
parametric dependence as before. The parametric estimation (4.16) thus remains intact. The number of the stretched 
degrees of freedom and the value of the gap need to be specified in (4.16) to describe the anomaly quantitatively. 

We start the quantitative analysis from discussion of the two-point instanton that is realized at 1 > 7 > only. It 
is an example of an instanton of the first type with both Se and A^; being of the order of unity (with respect to n). 
The dynamics of the two-point point instanton is characterized by the relative divergence of the points (drops) along 
with a simultaneous contraction of the sizes of the drops (see Section IIIA and Appendix A3). This means there are 
three different types of fluctuation degrees of freedom in this case. First of all they are longitudinal fluctuations of the 
stretched degree of freedom; second, fluctuations of the points (whole drops) in the d — 1 directions transversal to the 
stretched one; and, third, fluctuations of all the rest (2n — l)d contracted degrees of freedom (intrinsic fluctuations 
of the drops). One can calculate relative fluctuations of the drops and inner fluctuations of the drops themselves 
independently (it is easy to check afterwards that nondiagonal terms are negligible). Gaussian integrations account 
for relative longitudinal fluctuations of the points forms 



r/L] 



[l/2+AJ,']7/[2-7] 



2-7/22. 

DL7 



-l/[2-7l 



(4.17) 
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with as calculated in Appendix Bl. One finds that the d — 1 transversal fluctuations cannot be considered as 
Gaussian ones (attempts to restrict their study by a Gaussian level leads to divergence, see Appendix Bl). Hopefully, 
one can calculate the transversal Non-gaussian fluctuations explicitly. First, accounting for the relative fluctuations of 
the points (drops), at the initial (t = 0"*") and final geometries fixed, is performed by the method described in Appendix 
C. Second, one can account for the transversal fluctuations of and explicitly too, calculating a variety of rotating 
classical trajectories with T, R and r' taken from the rotationless trajectory, whereas the angular momentum ( 3.12| ) 



is nonzero (the trajectories are found from the auxiliary classical problem (3.8,^^9) but not from the full one (3.1,^.2) 



with a flxed form of the source function corresponding to the rotationless configuration). As a result, the trajectories 
with nonzero angular momentums give the same value of the classical action in the leading order in (r'/L)^, as for the 
directly stretched rotationless case (see explicit calculation for d = 2 in Appendixes A). In brief, accounting for the 
contribution of strongly Non-Gaussian transversal fluctuations, results in the r-independent multiplier (volume of the 
angular group). We discuss fluctuations of the drops themselves (the rest (2n — l)d fluctuating degrees of freedom) 
in Appendix B2. The fluctuations are sho rt-c orrelated, that results in r-independence of the respective contribution 
z/' into Z-^K However, the contribution (B7) shows an essential dependence on both T and n. Making substitution 



of ( 4.17 ), ( p7|) (Z-'' — Z^l^Z^), and ( 3.16| ) into (4^), and performing the integration over T in the saddle-point 



manner, one finally gets that characteristic value of the integral time is getting smaller, T ^ L'^ /[nD] ^ L'^ / D, with 
n oo. 

It should be stressed once again that the result derived for < 7 < 1 is not a consequence of a specially chosen 
initial geometry and source term - it is generic. For a majority of appropriate initial geometries and source functions, 
there is a special optimal configuration (that still may be difficult to find) of values and orientations of the initial 
momenta, making one distance diverge but all the rest converge dynamically. 

As it is shown in Section lllA and Appendix A3, at 7 > 1 there exist no alternatives to the common-type stretching 
of all the degrees of freedom: contraction of any distance (merging of any particles in a point) leads to a singularity, 
that is forbidden. And yet among all the generally stretched instantons it is preferable to get ones characterized by 
a stretching with at least one direction (dimension) kept stretching free (contracted): This is what we call "typical" 
instanton. "Typical" symmetrical instantons explained in Section lllA and Appendix A3 are 5*^ spherical instantons 
with m < d. Thus, let us apply the conducted above analysis for the case. The n-dependence of the potential term 
from the 2n particle Quantum mechanics ( 4.13| ) is estimated as Bij ~ 1/n^ for all the values of the particle index j. 



except of ones from a small vicinity (on the sphere) of i (every momentum is proportional to 1 /n in the dimensionless 
variables). For i and j being the nearest neighbors on the Sm sphere one gets [B]ij ^ ri7/[™-i]-2^ Qj^g (.^j^ drop all the 
terms beside the nearest neighbors if 7/[m — 1] > 1. Vice versa, if 7/[m — 1] < 1 one can replace all the matrix elements 
by ~ 1/n^ terms. All the kinematic matrix elements ICij are n-independent (strictly speaking for i and j being close 
to each other the matrix elements are even getting smaller with n ^ 00 than a constant, ^ n^I^ r]/[™-i]) Hence at 
7/ [to — 1] < 1 the energy characteristics are estimated as 5b ~ n^^, and A^; as a constant respectively. In the opposite 
case 7/[to— 1] > 1 (that is realized only if to = 2,7 > 1) one gets 6e ~ n^/I^^™"^)!"^, A^; ^ 7i7/[2(m-i)]^ Calculation 
of the n-dependence of the A term does not change the principal dependence on n of the energy characteristics. There 
is the extra parameter I which enters the anomalous answer and counts the number of typically stretched degrees 
of freedom. For the Sm instanton one gets I = 2(to — l)n. To conclude, contribution of fluctuations about the Sm 
instantons, d - 1 > to > 2, is estimated as - [r/L](AE+(m-i)n)7/[2-7][„2-7/22./[^^7]]-2(m-i)n/[2-7]. By this means 
the Gaussian correction appears to be of the same order as (even larger than) the saddle-point value ( |3.17 ), rendering 



the Sm saddle-points to be smooth out by the Gaussian fluctuations. Particularly, at 7 > 1 the contribution of the 
Sm instanton (saddle-point -I- Gaussian fluctuations) into the (e") correlator is getting smaller with Pe increase. The 
contribution is of no interest since it is negligible in comparison with the forced term contribution (possessing the 
normal scaling) that was dropped in the saddle-point approach from the very beginning. One recognize that the 
saddle-point calculus is not an appropriate tool for calculations of the anomalous exponent at 7 > 1. 

It was thus shown in the present Section that at < 7 < 1 the contribution of Gaussian fluctuations into (e") 
is algebraic (scale- invariant) with respect to the Peclet number (Pe = L/ra) and it is small in comparison with the 



"classical" value (3.17). The scaling exponent (2n of the scalar's structural function shows a flnite limit at n ^ 00. 

7 d + 2 - 7 



The exponent is calculated explicitly 
at < 7 < 1, Coo = 



2(2-7) 2(2-7) 



2-37/2+ a/16- I67+ 1772/4 . (4.18) 



There are relative fluctuations of the points in the two-point geometry that are responsible for the answer (4.18). 
At 7 > 1 all the saddle-point solutions discussed in Section 111 are smoothed out by the Gaussian fluctuations: The 
instanton calculus does not work in this case. 
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5. CONCLUSION 



It was stated in the Introduction that the idea of the saddle-point calculus is to make n the largest number in the 
problem. However, to establish the criterion of applicability of the saddle-point approximation at < 7 < 1 explicitly 
one should estimate contributions of Non-Gaussian fluctuations about the instanton and to compare them with the 
already found Gaussian corrections. If following the general scheme ( |4.4[ - [l.l4| ) to keep a nonlinear (say, third order 
over f term) one arrives at an extra factor y/T /[Rp{t)]''^'^ ^ Pe''^'^/^/n behind the dimensionless ~ s'^ term in the 
nonlinear variant of ( 4.13| ). The factor (along with the integral time T) is getting smaller with 71 —> cxd (the smallness 
makes the saddle-point become instant, and the "classical" action become, respectively, large). The observation is 
generic: all the higher order corrections to the Quantum mechanics describing the Gaussian fluctuations are small if 
n is the largest number in the problem 

>Pe'^/2d. (5.1) 



n 



It is remarkable that the anomalous scale-invariant result (4.18) has a more wider criterion of applicability than 
the method used for its derivation. Indeed, it is the zero mode of the primary operator C that makes the major 
contribution into the 2n-th order structural function. Choosing the scale-invariance to be the key for classification 
of zero modes, one finds that only n, d and 7 (but not Pe) can be entered into the zero modes. Thus, the anomaly 
(4.18) is valid even if Pe'^/^ is larger than n, however n remains much larger than d. It is important to note that 
(4.18) matches parametrically correct the perturbative results of the 1/d expansion at n being of the order of d. 

There existed two different expectations for (2n at n being sufficiently large. Decoupling of the molecular-diffusion 
term in the equation (that is not closed originally) for the structural function of the 2n-th order gives C,2n \/2ndQ2 
Another prediction is that (^2n tends to a n-independent constant determined somehow by d and C2 |p^,p|,p3[ . 
Thus, the recent work [ p3[ is based on an extension onto the passive scalar case of the method of operator product 
expansion suggested recently in the context of Burgers' turbulence |^^. Our calculations contradict the first prediction 
and support the prediction for <^2n to approach a constant at n — s- 00 if < C2 < 1- It is worth noting that an 
extended prediction of [ p3[ | gives us even more than an asymptotic constant behavior for the exponent at large n: 
C,2n was predicted to be a constant for all the numbers n having been larger than some rig. We cannot exclude or 
confirm the extended prediction here. It will require accounting for corrections to the saddle-point solution of the 
next non-Gaussian order explicitly. Let us stress that the equivalent problem in the Burgers' turbulence is not yet 
solved too. 

A comparison of the first-quantized formalism ( "Quantum" particles) presented in the paper with a second-quantized 
one ("Quantum" fields), that hopefully will be developed in accordance with the general scheme [Q, would be very 
instructive. 

It was argued phenomenologically 25 1 and confirmed quantitatively at d 2> n by means of direct expansion over 
finite velocity's correlation time [Q, that the scalar's exponents C2n are nonuniversal, they are crucially aware of 
the velocity's temporal characteristics. In the present paper we developed a theory for a peculiarly adopted (to the 
analytical study) case of the 5-correlated Gaussian velocity field. Nevertheless, the starting technical idea of the paper 
to replace in the operator representation by momentum of i-th particle in the path integral does not require any 
temporal or statistical restrictions on the velocity field being imposed. Hence it follows that it would be interesting to 
generalize the theory of finding C,2n at n ^ d for the more realistic case of finite temporal correlations and generally 
for arbitrary degrees of non-Gaussianity of the velocity field. It remains to be seen whether the scaling of the largest 
moments is nonuniversal in the limit. It is hoped that quantitative comparison of the future theory with experimental 
data (say, with measurements of temperature structural functions, up to order 12 [ p7[ ; for review of experiments see 
p8[) will be real some day. 

The theory prediction does not contradict simulations for rf = 2, C2 = 1/2 reported in Q (there, the velocity field 
was swept rapidly through the scalar to mimic the short-correlated feature of the Kraichnan's model). The largest 
tenth moment measured in the simulation gives Cio ~ 1.6085, that is smaller (as it should be due to the convexity 
inequality) than our asymptotic result Coo ~ 5.1. 
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APPENDIX A: SOURCE-FREE SYMMETRICAL INSTANTONS 



1. Si sphere geometry in d dimensions 



Consider 2n points equidistantly distributed on the I dimensional sphere: for example in the two-dimensional case, 
d = I = 2, we can use the following polar representation pk = (p, ^Pk), fk = {k — l) 7r/n, fc = 1, • • • , 2n, for the points' 
displacements. Let us consider the following setup 1) the angular momentum (3.12) to be equal to zero; 2) the initial 
[t = 0"*") spherical geometry to be preserved dynamically. Then, all the vector objects defined for some particle i, like 
Pj or Pi, are parallel to pi. The system of equations (3^ 3^) is reduced to 



1 



02 



Dp 



i^(2sin[(pfc/2]f ^((2d-7)sin2[(^fe/2]-d-l + 7) > 0, 



1 



ai 



2a2 



Dp2-i-^^ 2-7 

Thus, we arrive at the simple equation with the following boundary conditions imposed 

P 



d 
dt 

P(0) 



= r', p{T) 



a2Dp^-^p' 
{aiDp^-r + 2Kf 



= 0, 



R. 



Solution of the equation (A5) fixed by the conditions ( |A6| ) is 



Pit) 



dx 



y/aiDx'^-'y + 2k T J,., ^aiDx'^-'y + 2k 



dx 



(Al) 
(A2) 



(A3) 
(A4) 

(A5) 
(A6) 

(A7) 



in accordance with the energy conservation law ( |3.11| ) (to get the answer (A7) one could replace one of the basic 
equations (A1-A2) by the conservation law (3.11)). On the present instanton solution (A7) the action (2.11) gets the 
following form 



dx 



VaiDx^-'f + 2k 



(A8) 



The momentum of the k ~ O'th particle at zero moment of time is 

f 

P 



l^dx[aiDx^-^ + 2k]-^/^ 



(A9) 



t^o TVaiDr''^-"' + 2k 

Calculating ai in the two-dimensional case of S2 geometry, one finds that ai — [(4 — 7)^2/4 — (3 — 7)6i]/(2 — 7) and 

24-7r[(3-7)/2] 



r[(4-7)/2] 

2'^-^r[(5-7)/2] 
r[(6-7)/2] ' 



(AlO) 
(All) 



61 =^(2sin[^fc/2])2-^, 61 

k 

62 = ^(2sin[(^fc/2])*-^ 62 

k 

As n goes to 00 (and 7 bein g not too closed to 2^) ai goes to zero as n"'^^: If to replace 61^2 in the definition of ai 
by their asymptotic values ( A10| , All ), a remarkable cancelation, that gets rid of the linear over n term in oi, takes 
place. The cancelation occurs as a direct consequence of the incompressibility of the flow (the divergeless of the /C 
matrix). Thus, the major contribution into ai is estimated by the first term of the series over n. 
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ai 



(A12) 



where it is calculated that the angular size of the elementary cell, appeared after sectioning of the Sd sphere on n 
parts, is proportional to Tr/n^/^'^~^\ Considering ^2 geometry in higher dimensions d > 2, and generally Si geometry 
in d > /, one finds a linear growth of oi as n goes to oo. 



2. S2 geometry. General two-dimensional case. 



Let us fix the initial and final vectors, that define the symmetric geometry (say, pi(0) and pi(T)), to be nonparallel 
to each other. The situation starts to be more complicated, if one is going to solve the same instanton equation ([3.l|), 
as new parameters describing rotation of the vector come into the game. We consider here the two-dimensional case, 
that is the simplest one (rotation is Abclian in this case). In d = 2 the vectors pi and can be parameterized as 
pf — > pexp[i{ip + ipk)], pf — > pexp[i{'ip + ip + ^k)]i where scalars p{t), p{t) and angles (p{t), 1^(1) are time-dependent. 
Four equations describing the dynamical behavior are 



p+Dp p 



(3-7Cos2[i^])6i-t-l-i62 



cos[V'] 



= 0, 



iJj + ip+Dpp^ 
Dpp'^-f 



(7Cos2[^] - 1)61 + ^—^62 



sin[?/j] 



= 0, 



P = 



2-7 
Dpp^ 

2-7 



(3 - 7)&i + 



cos[V'] — 2 COS 



7 — 4 
bi + ^^b2 



sinhp] ~ 2k — sinhp], 
P 



(A13) 
(A14) 
(A15) 
(A16) 



where bi and 62 were introduced in ( A10| , All). The equations of motion ( A13| - AI6) are compatible with the con- 
servation law s of energ y (3.11) and the angular momentum (3.12) (that is pseudoscalar in this case). The system 



of equations (A13-A16) can be analyzed in full glory. Nevertheless starting from the point one drops diffusion term 



responsible for the ultraviolet regularization only. Then in the diffusion- free case, the system of equations (A13 
AI6 ) being rewritten in terms of the auxiliary x^y variables, {x,y) = /3'''/2(cos[7(y9/(2y^)], sin[7(^/(2y^)]), describes ; 
uniform motion of a particle with constant speed in the x — y plane: 



xit) = r"^/2 + A J^cos[7(^*/(2V^)]y^i?7 - 2i?'r/2r/7/2 cos[7(p,/(27^)] + r'^ 
y{t) ^^sm[-f^J{2^)]^W - 2i?^/2r'7/2cos[7V7,/(2V^)] +r'T, 



p'7/2 



where 



V = 



[(4- 7)62/4 
[(4-7)62/4- (3-7)61 



(A17) 
(A18) 

(A19) 



and the following initial and final conditions for p, ip dynamical fields are imposed: p(0) = r' , p{T) ~ R and (^(0) = 0, 
^p[T) = ip^. The solution ( [aT^ , |A18D means, particularly, that there are no equivalent trajectories for taken from 
the interval [(^c, | < 2t:^/^. Thus, considering the final values of and (p* -t- 2?™ from the interval to be equivalent, 
one observes y^/7-fold degeneracy (all those trajectories from the degenerate set are differed by the values of energy 
E and angular momentum M). 

The action S on the "classical" trajectory gets the following form 



4(2 -7)n 



_ 2r'T/2i?7/2 cos[7(/3 J(2^)] + r'7 



7'[(4 -7)^2/4- (3- 7)6i]i? 



T 



(A20) 



The momentum of the k — 0-th particle at zero moment of time (the object entered the self-consistency condition 
( |3.6|) and the pre-exponent term) has the following dependence on the initial and final conditions imposed: pi(t = 

0) =poexp[i(/?(0) + (pk +^"0], 
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7[(4-7)&2/4-(3-7)&i]^r 
X (^Wi-2r'i/^W/^ cos[7V3*/(2V^)]+r''' cos[7(^,/(2V^)] -r'^/^^ r"'/^-^ 

^ ~7V^[(4-7)6!/4'-(3-7)6i ]^ 
X y^i?7-2r'T/2i?7/2cos[7(^*/(27^)]+r'-r sm[-iLpJ {2^)]r'^' '^^^ . 



(A21) 



(A22) 



3. Two-point geometry. < 7 < 1 



The saddle-point system of equations (3T^,3.9) has some reduction feature at < 7 < 1: if we merge a group of 
particles in a point and choose the momenta (equal for the particles pasted together) to be parallel to the vector 
connecting the points at the initial moment of time, the problem gets rid of those superfluous degrees of freedom at 
all the latest times too - the group can be replaced by one particle. We will check the general observation for the 
two-point case in the present Appendix. 

Here we consider a geometry formed by two groups (labeled by -I- and — ) of particles (n+ + ri- = 2n) each merged 
in a point with a position p+{t) or P-(t) respectively. Indeed, there exists a solution of the saddle point equations 
( |3.8| , |3.9| ) that preserves the symmetry dynamically. That is specific about 7 < 1, it is an algebraic decay of JC"^'"^ 
when r goes downscales; one can put the particles in a point without any divergences appearing. The number of 
dynamical degrees of freedom is reduced from 4nd to 2d. One gets 



pI + 2Kn±pl - n^K:fv% = 0. 



(A23) 
(A24) 



Coming from P-|-,p± variables to collective ones p = ?i+p+ = — n_p_, p (the full momentum is constant in the 
reference frame associated with the center of inertia) , one gets the system of equations that does not depend on the 
numbers of -|- and — particles at all 



V 



4kp" -I- 2/C^'3-'' 



(A25) 
(A26) 



Notice that the rotationless variant of (A25 A26D is transformed to (Al,A2) and the two-dimensional variant of 



2(d-l)/(2-7), and k 2k and t] 7^/4 



( A25| , A26) is transformed to (Al2-Alfi), if one performs a reduction ai 
there. Thus, the formulas for the "classical" action (A_8) and the initial momentum (A£) are valid with appropriate 
changes (the multiplier [2n]~^ should be accounted addit ionally in the action) in the case of the two-point geometry 
too. One finds that the effective "classical" action ( ^.17 ) is n-independent in this case, if the values of T and R are 
considered to be fixed. 

To check the dynamical stability of the two-point geometry let us consider a "dumb-bell" geometry (floated variant 
of the two-point one) with the characteristic size of the drops of the "dumb-bell" being initially (at i = 0"*") much 
smaller than the relative distance between them. The major question to ask is: Will the "dum b-bell" geometry be 
preserved dynamically? To answer the question, let us go back to the classical equation (^.9|). Introduce a small 



excursion of a particle ps from the drop's center. Then, keeping the leading (over Sp/p) term in the rhs of (3.9) one 
arrives at the following equation for 6p 



Sp°' = pp^'-f {Sp°' - d{Spn)n°') . 
Substituting the already known law of the basic two-point stretching, one gets, 

(2-7)/[2(<i-l)] 



5p±iO+) 
Spwit) 
6p\\iO+) Vr 



r 

P. 

p\ 1-7/2 



(A27) 

(A28) 
(A29) 
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where 5p±^\\ are transversal and longitudinal (with respect to the direction of the "classical" stretching) sizes of the 
drops. To conclude, the ratios of the drops' sizes to the distances between the drops are getting smaller with time 
and the two-point geometry is indeed preserved dynamically. Note that (A28) is a "classical" manifestation of the law 
of volume conservation valid at 7 = (the respective law of area conservation allows solving the problem at d = 2 
explicitly |l3|Jl9| ]): p[Sp±]'^~^ does not depend on time. In other words, (A28) is a statistical descendant of the volume 
conservation law (incompressibility condition) valid for any particular flow. 

Considering the dumb-bell geometry at 7 > 1, one finds that it is destroyed dynamically: The size of the drops and 
separation between them, in the initially served two-point geometry, turns out to be of the same order at the latest 
times. 



APPENDIX B: GAUSSIAN FLUCTUATIONS ABOUT THE TWO-POINT INSTANTON 

1. Relative fluctuations of the points 



Consider Quantum mechanics (4.11-4.14), that appears when accounting for the relative fluctuations of the drops 
only. Thereare d essential distances in the case: f = r^iy — f{2}- Operator C' ( 4.13| ), rewritten in terms of s (related 
to f via (ITTol)), IS 



d+l-7„^ 



7 



4-7,2(2-7)^^,„^2(2-7)^^^^^^„ 



(2-7)^ 



+ 2-7/2 + 



d-1 
1 

0?+ 1 -7 



V"- 



(ns)^ 



(Bl) 



where n = p'^'-fp'^K It is seen clearly from (Bl) that the transversal part of the effective potential is not bound from 
below. Accounting for nonlinear terms is required to regularize the divergence (see Section IV and Appendix C for 
an explanation of how to avoid an explicit calculation of the terms). By this means, to describe the longitudinal 
fluctuation one should make s be parallel to n in (Bl) and so deal with the reduced operator 



d-1 d^ 
2-7ds2 



{d + 2--f){d-l)l d 4-37 d (3-7)(2-7)2 
2-7 'sds ^ 27 ^ds 2-f^{d- 1) 



(B2) 



In the long time (t ^ 1) asymptotic the major contribution into '^t stems from the lowest eigenstate of C" . The 
eigenfunction of the ground state is an exponential of the quadratic form 



gr 



exp 



as 



(B3) 



In the Heisenberg representation the eigenfunction appears as exp[A£;T]^'gr, where the ground state energy for the 
eigenvalue problem gets 



Ab = 



27 



2 - -7 + v/(2 - 37/2)2 + 2(2 -7)(3- 7) 



(B4) 



2. Inner fluctuations of the points 



/C, A and B responsible for the inner fluctuations of a point (disk) are estimated as 

/Ci-V"^ ^5~pV-^ Bs^p'Sp-^ 



(B5) 



where 6p{t) is related to p(t) via (A28 A29). Thus, instead of (4.10) one has to perform the following transformation 
to the dimensionless , s^^ variables 



7 ^ 
2X' P \ p 



1-7 



Sssi = 



^jT5pp^~^/^ 



(B6) 
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to get rid of the temporal dependence of the Hamiltonian. The resulting Quantum mechanics will yield positive 
spectrum with a constant gap. The dependences of A^; and Se on n are estimated as const and ~ 1/n respectively 
(the potential /C-term has an extra smallness ~ ^/n^, while the number of elementary excitations is ~ n). The 
dimensionless time of evolution is short ^ [i5p(0''")/p(0"'')]^~'''/[l — 7] ^ 1, if 7 is not too close to 1. These collected 
observations result in the following asymptotic for ^f^r (analog of ( 1.15| )),^^t Sx (f [L/r']-'^"'''/^) (we have calculated 
( A29| ) here). The respective contribution of the fluctuations is thus scale (r) independent 



n 

drop 



exp 



dct[g]\ 



Sx{i[L/r' 



'11-7/2 



2dn 



q = ma.x[l;L'^/[nDT]], 



(B7) 



where, accounting for scale dependence only, we use Sx{f[L/r']^ ''^^), as a function of f, decays on the scale 
^/i-7/2_^7/2 _ rL'^/[nDT]. 



APPENDIX C: INTEGRATION OVER THE SOFT ROTATION MODE 



The measure of functional integration in (2.1C) is quasi-invariant with respect to a slight rotation of the fields 

p,{t)^U-'{t)p,{t), p,{t)^U-\t)p,{t), det[J7(t)] = 1, J7(T) = ^7(0) = i, (CI) 
where U{t) is a unitary matrix realizing the d x d rep resentation of SU{d). Quasi-invariance means that an extra 



term appears at the transformation of the action ( 2.11 







AS= [ 




Jo 





a/3 



P?P^dt. 



(C2) 



To fix the gauge quasi-invariance we will use a general method p9| ], that is popularly known in field theory. To 
integrate over the soft mode we will put under the functional integration ( ^.10 ) the unity 



1 



VU{t)S 



{iipi)' X iip^ 



p{ det 



'su 



(C3) 



where, ni = pi/pi and x stands for the antisymmetric vector product of d dimensional vectors. Thus, we use a 
requirement for one of the particles (labeled by 1), to not rotate over the origin, as a gauge condition (without loss 
of generality one can choose an arbitrary direction, chara cter ized the trajectory, to be nonrotating) . For the sake of 
simplicity, let us consider the two dimensional version of (|C3| ) 



1 = J Vip{t)S 
U{t) = 

Direct calculation gives 



X Upi 



pi det 



(upi)' xUpi 



pI 



Sip 



(C4) 



cos Lp sm Lp 
- sin Lp cos 



(Upi^ xUpi 



Pi 



+ [Pi X Pi]/pI. 



(C5) 



Thus, the determinant from the rhs of (C4) does not depend on the dynamical field pi, and we can drop the 
determinant. As a next step let us perform the change of variables (CI) in the original functional integral (2.11) (with 



the unity (C4) substituted into the integrand). The Jacobian of such a transformation is unity. The S- function from 



the rhs of (C4), describing the gauge condition, turns out to be (yS-independent. We get finally that the only factor 



calculated integration over the soft (gauge) mode is 



yp^exphA^M], A5 = y ^e^Pptp^dt, 



(C6) 



where i is the antisymmetric 2x2 tensor. Further, the integration over Vip, being performed, reduces (C6) to the 
(5-function 
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(C7) 



The condition under the (S-function, which is satisfied on a saddle-point solution as describing conservation of the 
angular momentum (3.12), is thus valid for fluctuations too. 
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